y 


AD-A04% 674 FOREIGN TECHNOLOGY DIV WRIGHT-PATTERSON AFB OHIO F/G 9/3 
PROCEDURE AND PROGRAM OF THE CALCULATIONS ON TSVM (UBM = DIGITA=-ETC(U) 
NOV 76 VV ILINICHNINe L A KOROLEVe L S LEVIN 

UNCLASSIFIED FTD=ID(RS) T-1606-76 NL 


FOREIGN TECHNOLOGY DIVISION 


PROCEDURE AND THF PROGRAM OF 


PROCESSES IN 
CIRCU 


AD-AOF4¢ 674 


FID- 


i 
Se eee — 


FTD-1D Nove 


PROCEDURE 
TSVM [UBM - DIGITAL COMPUTER] 
PROCESSES IN LINEAR ELECTRICAL CIRC 
CONCENTRATED PARAMETERS 


AND THE PROGRAM OF 


CA 


Levin 


59 


English pages: iz 


Issledovatel'skiy Institut 
Blektroenergetiki, Izd-vo 
Moscow, Nr 395 1971, PP. 3 


Country of origin: USSR 
This document is a machine aided translation. 
Requester: FTD/ETDP 
Approved for public release; distribution unlimited. 


THIS TRANSLATION IS A RENDITION OF THE ORIGI- 
NAL FOREIGN TEXT WITHOUT ANY ANALYTICAL OR 
EDITORIAL COMMENT. STATEMENTS OR THEORIES 
ADVOCATEDOR IMPLIED ARE THOSE OF THE SOURCE 


AND DO NOT NECESSARIL Y REFLECT THE POSITION 
OR OPINION OF THE FOREIGN TECHNOLOGY DI- 
VISION, 


FTD- ID( RS) T=-1606-76 


TR 


BY: Veo Va Dbl rtehncin aA emOmol ewes. Ws) 1S. 


Source: Trudy Vniie Vsesoyugnyy Nauchno- 


UNEDITED MACHINE TRANSLATION 


mber 1976 


LCULATIONS ON 
ANSIENT 
UITS WITH THE 


"Energiya, 
Frame lies < 


PREPARED BY: 


TRANSLATION DIVISION 
FOREIGN TECHNOLOGY DIVISION 
WP-AFB, OHIO. 


Date 17 Nov jg 76 


U. S. BOARD ON GEOGRAPHIC NAMES TRANSLITERATION SYSTEM 


Block Italic Transliteration Block Italic Transliteration 


Aa a A, a Piece Pp Ry of 
67 26 6 Bin 19, { ¢ aR 
B 8 é V, Vv T fee 


m 


A P 

5 Cc 

B T 

Is y 

q ) 

E x 
HOOK KK Zi, 2h ae LU Ayes as! | 
33 Sos Ley Bi heel bef Ch, ch | 
uu Hou Bt wow Sh, sh | 
Ww Aoa any Uy uw Ul wy Sheh, shch | 
Hoo K x Ke & Bt a 

fla Weve > bl ow Yee 

M Ly am b b bb ' 
HH Hw N, on 33 3» E, e 
Do O o Ds 30 1 fo » YU -¥ 
on i tk Baaw AR eS Ya, ya 


*ve initially, after vowels, and after s, b; e elsewhere. 
When written as 6 in Russian, transliterate as yé or eé. 
The use of diacritical marks is preferred, but such marks 


may be omitted when expediency dictates. 


Alpha Ay es @ Nu , y) 
Beta I 8 KL 

ramma Y micron QO oO 
Delta A iL | c | 
Epsilon BE € Rh ge 40 e 
Zeta G C Lema } 0 S | 
Eta ih) Tau ot 
Theta 60 6 §$ Upsilon YT wv | 
Tota [ 1 ato! >’ @ 6 

appa K nN F x Cni \ x 
Lambda A \ si . wy 
Mu M ou Omega RQ w 


ud 
i 


'D=ID( RS )T=1606-76 


RUSSIAN AND ENGLISH TRIGONOMETRIC FUNCTIONS 


14 . 
English 


NUSSLIaN gf 


Q 


cacihi 


GRAPHICS DISCLAIMER 


All figures, graphics, tables, equations, etc. 
merged into this translation were extracted 
from the best quality copy available. 


19-12-7/ 


MEAS 


“SUB: 


the 


LE VNEEE [Haine = 


‘OWER ENGINEERING]. 


hyn, ie er 


i transient 


>4 77 


ALL— UNI 


Al1-U 


lectro-energeticists. 
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Lansient processes in li 
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‘ Ww 
PARAMET eet. 
design of electrical s 
tudy of electromaqnetic 
lse in 


Operational pract 


x9’ damages cn power tran 


nj and the 


selection ot 
nmrocedure and were carri 


processes 1n 


GON SCPENTLIEIC RESEARCH 


nicn scientific research 


lations on TSVM f[y&™" - digital 


near electrical circuits with 


ystems frequently appears the 
translent frocesses. A Similar 
ice for the analysis of the 


snissicns of high and ultrahigh 


equipment, etc. In VNIIF was 
ed out the program of the 


digital 


cemputers (TsSVM) in 


LL eee ys . 


distr} ted circuits (concentrated linear reactive cell/elements are 


repla the cuts of long lines) [l. 1]. This work headed by 
docto: fech. ScLences, by Prof. A. I. Doloincy, wader 
Bana) ¢ innal of whom beqan to be developed cther directions of 
the c4 itlitions of electromagnetic transient processes on TSVM. one 
of su irections was use TSVM for the calculations of transient 
Cur re ind voltage/stresses in electrical circuits by the 
concen? “=a cellfelements (cell/felements ct grid/network with the 
GSS tri rut parrmeters are represented by the appropriate replacement 
Sepen 

‘nt atticle are represented the frocedure and the program 
Of “tl leulations of electromagnetic transient processes on TSVM in 


electrical circuits with lumped parameters with the use of numerical] 


methcds of integraticn. 


Fubdamenta]l principles. For the calculaticn of electromagnetic 


Fans} oht roceSsesS on TSVM by the numerical methods of approximative 


tantegqrat must be comprised the system of differential equations in 

can¢ EM. | 
| 
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Fort ntasition cf the reference system cf differential equations 

We uti ; equations of the junction/unit potentials, which 

5 ' 
efarac*-1 this electrical circuit: 


ee Oa PAGE 
muh, (1) 
where 
a | matrix cclumn k cf the potentials of the assemblies 
Uz 
wale. 


circuit with respect to the frase line assembly, 


ay Jtential of which is accepted equal to zero; 
fa || ; : foree 2 
rere ae matrix column k of the assigning currents in 
= 3 


lll the assemblies cf replacement scheme; 


1 gucgas s+ —ga || the square matrix/die of order k 
— Ul Cae = s 
e=. But on — 8s | its own and vzamnykh conductiv 25 
Gar Sse > + Bow ll the assemblies of circuit. 
It we consider that 
t 
d ° 
P="! piss 
0 
that rc tor/resistance of kranch “an between two assemblies 
and n i xressed by the following rating value from p: 
4 A 1 = —l © 
Lain P) = PLitin + Ring + PCT (2) 
wher: L—, R= ,Ca- are values of inductance, of effective re 
mn mn mn 
and inv capacitance/capacity of the cellyelements, connected 
seri: t issemblies m and ne 


Bmn (P) matrixydie g (p), in line nm 


mM 


istance 


in 


Se be? PAGE 4 
cenfrortiny and column n, in oShchem form is represented by the 
expre t 
Sper 

U = EUnm= eee ee eee gh 

Gane ee EP Dhigat Seah” Oma 
in which ¢ amount of terms it depends on the number of parallel 
conn between assemblies mand n (m 3 n)- and also the 
number connections of danogo assembly (m = n) with other k the 
assem! 1i xf the Snreual carcuat.. 

(1) the unknown potentials of assemblies are determined 
throu t known assigning currents i and matrix/die z = g-! their 
own rutual impedances of the assemblies of circuit, i.e., 

u(p) = Z(p) 1. (3) 

; ; : : : : | 
Matrix Z {p) is calculated by matrix inversicn g (p) and as it 
will ywnh below, it can be represented in the following form: 

ay. (p) Qia(p). » + ain (P) | 
1 Qa, (P) @aa(P)+ + + an (P) (4) 
CUP Spy | sa ee gon we 
ny (P) Ana (P) © + + Ann (P) Ht 
>in Z til matrix/die 7% (p) - the correct rational 
fract ; Wii Le and numerator Ann \P) and denominator A (p) - 
pel yt ) Ds 
. Ann (P) oF OP He gD? te (5) 
Zmn(P)=—S(p) Bet dip+ dpi +. + + 
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he 
rte 
re 


with 


asseul 
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CCrres 


matrix 


equal + 


that 


equafi 
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(3) ks by account (4) the volitageystress of any assembly o 
t to base line is equal to: 
ban PD (p) P . ‘ 
tm id a ene 3 (p) Las (6) 
n=\ 
larger interest are of not the potentials of assemblies 
t+ to base line, hut the potential difference bet ween 
m and j, which it is easy to chtain, replacing in (6) 
Amn (P) by the coefficients, equal tc a difference in the 
nding coefficients of the numerators of lines m and j of 
lie 2 (p) (4), Lees, 
Fn (P) = in (P)_ 
Um4 (P) ~Y) fanlel tine) Se a (p) “ inte (7) 
nazi 
“urrent in branch mj with resistcr/resistance Zmi (P) is 
atid (P) 8 
Imi (P) = lnl\P) v (8) 
+ termination of currents Ing (P) it is possible to use 
(7), after replacing 4 (p) with value &4* (p), equal to 
"(p) = 4 (p) 25) (P) 
a ) as ohh ) . 
Fatt py aYi taal mn aaa .. (p) b. (9) i 
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LG MUS DE oted tha equ ations C6), CT (9) 1 wri ¥ 
Single pattern; 
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=yr- 0 
y= )) Bip) ” 
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Ee 1 \ 1 ; : 
wher y are unknown value, and A(p) and B(p) are polynomials 
respectively of reer gq end VILO eoef fic rv T fey 
-alecul 1 fr known (oreset) parameters of nitial 
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[ec operating equation (10) corresponds the following differential 
equaticn: 


dy dry dtu 
Buy + Biot Bag tt Be ae 


k 
bein dal, 
= Sau ftag os + Ain Gt + Agn dia ava (11) 


n=l 
solution to differential equation (11) is reduced to the 
soluticn of the system of differential first-order equations, which it 
is not difficult to cbtain vvediyem the supplemental equations, which 


reduce the order of equation (11): 


Gaya 
duo 


a Hi 


dyy 
“de aye 


. } (12) 
dys 2 
(iy eb 
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ate ene + Ain age Amn ~qae bt Aan iam | 
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Figs. 1 To the circuit design, which contains the inductively 


connected cellyvelements. a) the group cf the ends of chain, connected 


induct ively; b) replacement scheme. 


Key? (1). Complete polygon. 


! Meal Ay prey 
ip I Noanoi ui . 
IS Ser 

J eee t 1 ' r Mn0eg0- 
, ‘ K YEOMOHUK k 

| ! Zan! PO 2galan! Zpphn! ! y i 
| n ) a. n oa ON 

a. 
64 


PAGE a ® 


solution cf obtained system (12) is conducted hy the 


methods of numerical intearaticne 


sootnovenlya make it possible tc calculate the transient 
current and voltage values in schemes at the absence of inductive 
coupl i petween cellfelements. The use of a wethod of junction/unit 
potentials in the presence of vzaimoinduktivnykh connections is 
possitle with the application/use of a decoupling using the method of 


cemel 


use this method for the calculaticns of transient 
proces lor us consider that the internal resistance of branch in 
the ini l scheme (Fig. 1A), connected Fetween assemblies j and j* is 


egua! 


a = pL. 
2, (p) = plz, 


spectively mutual impedance retween the vzaimoinduktivnymi 


Ss Kk* and 337 is equal tos 


Lae (Pp) = PMagy jjr- 


jetermination cf complete polygon {Fig. tb) first 


matrix/die Zss the inductively connected 
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u 
ther is located matrix/die Yass tc reciprocal matrix Z,;, 
Va pale 
w ieasignate matrix element Noy by ¥,, Where j, k= 1, 


»» by, then to the reactivity of the branches cf complete polygon 


| 
Lye 5 (13) 
| 
| 
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-rlculation of currents in the induct ive-connected 
nts of replacement scheme is carried out as follows: 
| —i 
ln=—p Zu yj (14) 
I is a matrix cclumn of the unkneowr currents; uy), - the 
RO jj 
ylumn of potential differences cn the end/leads of branches 
tering the group of vzaimoinduktivnykh cellyelements. 
itufing in (14) values u 


yo calculated on (7), we 


pression for the calculation cf currents in the 


}. nts, connected by mutual induction, analogous (10), 


n is form/shaped the system of differential equations. 


‘aleulation of matrix elements z (pf). For the calculations 
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consect 


Luact 
each s 
ket we: 


cases: 


elem: 


Pade 


ccelu 


]+ 


nad to this method complete 


xX @alements Z (p) 1s utilized by a logician of the method of 


shchivaiya, utilized for the calculaticns of the 


thesikh comprising short-circuit currents on TSVM [l. 3, 4}. 


network is formed by means of the 
ive addition of kEranches and calculation of matrix/dies 


for the being obtained sequence cf nastichnykh circuits. In 


qe of growth is included one branch sq replacement scheme 
issemblies s and gq, resistive ye Here they can be two 


he newly included tranch sq introduces the new assembly 


9§ partial matrixydie increases per 


UHDT) . 
newly included kranch sq Joes net introduce new 
(poyadok of partial matrix/die does not increase). 

respect to these two cases the calculation of matrix 


Zaact is conducted in different ways. 


first case the calculation of the cellfelements of new 


line 1S carried out on the following cf the formulma: 


a) m#4q, Ziq = Zan = Lom ae (15) 
6) m= 4, Lag = 2a +25. < 


special when 


branch Sq is connected to base line 


case, 
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alm #4 Zing = Zan = Zam =, (m= 1,2, Q—I)i } (16) 
6)m=4q, 2q=25. 


second case the calculation is ccnducted into two stages; 
-omputed the cellyvelements cf supplemental column and 


a) m#L, Z mt, = 21m = 2me— 2mg (n= 1,2, wo b— li (17) 
6) m=L, 2p, =2Zyt 2qq— Mea + 25. 


special case, when branc} sq is ccnnected to base line 
, Ld ~ = O, 
a)m#L, Zin, = — Zmg (m= 1,2, 00, L— I); ) (18) 
6) mal, 2), 2g t+2y, 


eliminated newly introduced column and line l, i.e., we 


unknown cell/elerments of the new sartial matrix/die: 


2 mt. Zin 
Zin = 2 mn sea gaa (19) 


number of line; n —- the number cf cclumn; L - the 
supplemental the lines also cf cclumn; 2. = ‘he 


of old partial matrix/die (withcut the account of branct 


ae is a cell/element of new partial matrix/die (taking 
branch 39). 


nto account (2) and (5) formula (175) = (18) they will be 
tollows: 
I, s #0 
Gama (P), 
a) m#q, Zmq(P) = Zqm(P) = Zma(P) = A(p) , 


Qua (P) + 4(p) 25 (P) . 
6) m==q, Zyq(P) = 4 (p) : 


(20) 


s= 0 


5h ia PAGE 


0 
a) m# 9. 2mg(P) = 2Zam(P) = TT py? 


A (p25, (P) . 
6) m=, Zaq(P) = cases 7-1 Uae 


(21) 


Page 4 Dig. 
HW S40 
a) m# Ly Zins (P)= Lym (p) = ee 
7 ) (22) 
Ags (P) 7A gg (P)— 2a (P)+4(p)Zz (Pp), 
6) m=L, Z,, (p) om ———— ’ | 
s=0 
ame Z mt (P) = Zz (9) = — SD; | 
Maq(p) + 4 (p) 255 (p) f eas 
6) m=L, Zi (P) = eo hae se | 


, the application/suse of formulas (20) - (23) is reduced to 


addition (vychitaiyu) and to the multiplication cf polynomials from p. 


basie difficulty in the use of a method of growth conformably 


to th ‘alculation of transient proyesscv consists of the 
recalcalation of the cellyelements of partial matrix/die Zuac. with 
axcept i limination newly introduced stol'tsa and lines 1 according 
te formula (19). At each space of its use cf poyadok of numeratcr and 
denominator of fractions (5) in limiting case it can increase 2 times. 
This it is possible to avoid, by converting the algorithm of the 


except ion/e limination of the newly introduced column and line L. 


Ye5ootVetsStvil.s (22) in the general case we haves 


a —_ ee 


Where 


acc cur * 


sq). 


has as 


where 


A°(p) 


re 


finally instead cf (19) we have: 


Ginn (P) 
Zinn (BP) = Riv (26) 


" ann (Pp) Au (P) — ahi (Pp) an (P) 
Ginn (p) = Lae MS es 


Gy, (2) + Aqq (P) — 2ayg (p) + 4° (p)Z55(P) — aH(p) 
i) ue ne el aerate 7 etaaaa mAs (p) 1 (24) 
1S a polynomial A (p) old partial sgatrix/die 
tranck sq): 
A" (p) = ays (2) + gq (P) — 2g (p) + 4° (2) 255 (P) 
1 A (p) new partial matrix/die (taking into acc 
acite (19) taking into account (24) and (5) 
z" (p) = ee (p) _ Ant (p)a’.,(P) pa 
mn Ac (p) ae (p) AN (p) 
ann (P) 4" (p) = inp (P) a5 n (P) 
ae “2 c u ie (25) 
Ac (p) A" (p) : 
rossible to show that the polynomial 
@(p) = Ann (P) A" (p) —a5,1 (P) af, (P) 
iiVidetC/denominator polynomial Ac (p). 
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Input routine and the initial data 


kKtivnye connecticns? (Ss Yess (4). 


(2). Are there 


NOe (5). 


f vzaimoinduktivnykh connections using ¢he method of 


One (6). Calculation of the matrix/die of its own and 


lances and its printing according ¢ 


t 


sf the system of differential eguaticns. (8). Solution 


of differential equations. (9). 


Frinting of the 
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iS 

. after the switching on of all branches of the circuit from 
of on: ' consecutive calcnlation of the cell/elements of partial 
matri> i finally we will cbtain matrix/dies Z (p) in form (4). 

rt characteristic of program for TsVM m-220. The initial 
data ¢ t | calculations of transient processes on TsVM are the 
Compor: vVilnes of the Substitutions, which include coefficients of 
Des t the operating fresistor/resistances of the branches of 
sehene, roptficients of mutual inductance of network elements and 
the Gc lcients, which characterize value and the form of the 
assi currents in assemblies. In this case the determination of 
the u values is carried out by the zero initial conditions and 
the p1 © Of the assigning currents under each assembly of the 
hie as tk ircuit, which take the following form: 

bm = Atm S10 (Wot + Ym) + Aime o™ ’ 

wher: Atm, Asma) ime Orn are constant cceffticients, 

‘ notation of these data are assigned special places after 
prelj ty coding the assemblies of circuit by the numbers of natural 
Ser i« fter the introduction into working storage TsV¥M of the 
indicat initial data and data on the character of the unknown 
infctr the calculations of transient precess occur automatically. 
Page 


veloped program (Table 1) cotisist 


Cc 


f three basic building 


coett 


fesist 


diffe: 


enconu 


VS (pre 


matriy 


decr: 


ind vc? 


+ 


+} 


e block/moduleyunit of the replacement of 


Lynykh network elements by equivalent replacement scheme; 


eo 


i: 


“le 


, 


sembly of the calculation of the values of the 


ot the matrix elements of its cwn and mutual operating 


Stances; 


sembly of composition and Solution of the system of the 
y i 


juations of the unknown values. 


bly of the replacement of vzaimoinduktivnykh network 


c 


signed for five groups of the cellyelements, which have 


7*-10n. Each such group can into its cchered™ contain five 


connected inductively». Since during the calculations 


nrocesses, especially in the secondary circuits, can be 


more complex forms cf inductive coupling, in proaram 


he possibility of an increase in the maximum order of the 


Zo, Which characterize each such group, because of a 


amount of groups of the cellyelements, connected it is 


jecoupling of vzaimoinduktivnykh connections the 


fer/feonverts to the construction of the matrix/sdie of the 


own and mutual impedances cf replacement scheme and the 


the coefficients with operatcr's different deqrees p. 


‘ighed for the schemes, which ccntain to 40 assemblies 
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and # ches with the maximum degree of tte pelynomials of 
Cationil-tractional functions (5), equal to twelve. On those which 
were o>} rlated elementammatritsy Z (p) and the initial information 
are for taped the necessary systems of differential equations for 
the det ersingtion of the unknown values, which are solved by runge - 
Kutta '* riod the Sth order. Is simultanecusly feasible the 
calculirion to 15 unknown values. 

cencl terivations. 


ievelopment of this program showed that the proposed 


broc: can be successfully used tor the calculations of 

electi ‘gnetic transient processes in ccmplex linear electrical 
circu: “ith lumped parameters on TSVM with the automation of all 
stag the calculation. 

r fundamental positions of the develcped algorithm make it 
poss) *o propagate it to electrical circuits with the large degree 
Gf dif tial eguation, circuit, which cortain nonlinear 
cell + 
Appl]: impendix. 


le of the calculaticn cf transient process. Let us 
eXd th) calculation of transient current during triphase korkotkom 


closi rting on line 500 kV. The netwcerk cf primary network is 


b. 


KV, 


nlbic 
INLIC 


nted ¢ 


by 


from 


Triphase short 


length 300 kno, 


two assemblies 


circuit 
which is 


500 kV, 


generators and transformer 


he following parameters cf 
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occurred at the end/lead of 
lccated under 

each cf which consists of 
group. In the calculation 


instrumentation. 


L=300Km 


Fide. 2a » the calculation cf transiert process during triphase short 
circui line 500) ky. a) the scheme cf grid/network; b) replacement 
sche 

fable, 2. 


Bassnvenas | QD Conportnsaenie verven 
/ peTuen | | 


as Pestle 
0! | 0 2 | 0 
12 i 6 0 90 
23 ) 3,24 45 
03 | 870 0 0 
34 ie 3,24 45 
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Key: , signation of branches. (2) le Fesistor/resistance of 


0.021 ) 
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Pours Wi =135 000 KVA, Uw=13800 in, x4 =0,135, 
ss Of the transformer group: Wu =135 000 KVA, 
525000 | : ‘ Meare 
Uae | 13-800 in, Uxn=135%. lep 500 kV: Ps = 300 km, Ip = 
, & = 043 T7Km, Yor=3,87 «19-8 S/Kt@. 
zt) depicts the calculated replacement scheme, in which 


(ge; - 


lies and fa - the earth 


Approximately in the 


- m=05 Ohms Value enf 


by identical, 


| 


equal © = 


calculation on TsvM, the 


3 eG The versicn of this 
py the selected 


Lbileaw 


esies + 
ances Of 


the equivalent ind 


culeulataon 2£t 


4€0°10, sine r (7 - 


numbering 


Eranches in 


and ohmic resiStance of 


resistance of the duct of 


was accepted: 
all 


generators was 


time in 


assemblies of replacement scheme 


numbering is shown in Pig. tbe 
in TSVM will be brought in 


the following form. 


0. 
aft, Sequsnoepa mma 


Pacvernnan xpubas 


Vasa feet ‘alcnlated and experimental curves of short-circuit current. 


Key: ‘ iscillogram. (2a Cabkcibated curves (3). 2S glad. 


PAGF 
lation for assent] 1 must be alse carried the 
auplira nd the initial phase of source cf current IT = 18 0103 Sine 
T (dine l8O XK 108, Gu=O), 


e and also the value ct angular frequency (Wo = 


ve Calculated curve ef short-circuit current at 
carried out by zere 
ectness of the obtained 
seccndary current of 


eStablishedsirsta : > this line. 


Fa 


cte the satisfactory ajreement 


9stsillcgrammnoy curves, 
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